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Abstract 

We propose a novel molecular device that pumps heat against a thermal gradient. The system 
consists of a molecular element connecting two thermal reservoirs that are characterized by different 
spectral properties. The pumping action is achieved by applying an external force that periodically 
modulates molecular levels. This modulation affects periodic oscillations of the internal tempera- 
ture of the molecule and the strength of its coupling to each reservoir resulting in a net heat flow 
in the desired direction. The heat flow is examined in the slow and fast modulation limits and for 
different modulation waveforms, thus making it possible to optimize the device performance. 

PACS numbers: 63.22.+m, 44.10. +i, 05.60.-k, 66.70.+f 
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I. INTRODUCTION 



A heat pump is a device that transfers heat from a low to a high temperature reservoir 
by applying an external work that modulates the system's parameters. This paper discusses 
a molecular machine of this kind. The analogous electrical device that transfers charge (or 
spin) against the electrochemical potential bias was studied theoretically 1] and demon- 
strated experimentally in an open quantum dot when varying both the dot voltage and the 
tunneling barrier heights Q|. 

In a prototype particle pumping machine, each cycle begins with isolating the system from 
one reservoir by reducing its coupling to the system, while applying a potential that drives 
carriers from the other reservoir into the system. Next, this configuration is reversed, the 
system is coupled to the previously disconnected reservoir and isolated from the previously 
connected one, and its potential changes so as to drive carriers from the system into the 
connected reservoir. Consequently a net current is flowing through the system. A basic 
requirement for demonstrating pumping operation is the modulation of at least two internal 
)arameters. Out of phase modulation provides an adiabatic (reversible) pumping operation 
jij], while in the general case, quasiadiabatic (irreversible) processes can be realized 

Motivated by the growing interest in nanomechanics and quantum thermodynamics (f| , 
we present here a molecular model for a thermal pump that is based on similar operating 
principles. Other thermal devices that have been envisioned recently are a heat rectifier 
M a tana, tr an Slstor , Q aad even a m een ml an a log „ t a laser Q. A S with my 
machine, one seeks optimization of performance with respect to both efficiency and power. 

In our model a molecular unit connects two spatially separated left (L) and right (R) heat 
baths held at different temperatures, and transfers heat from the cold (C) (henceforth re- 
ferred to as the left side) into the hot (H) (right side) reservoir. An external force modulates 
the energy level structure of the conducting molecule and consequently its effective coupling 
to the reservoirs (thus providing a modulation of two system parameters while modulating a 
single physical variable). This system is shown to operate as a heat pump that can transfer 
energy from a cold to a hot reservoir. 

Similar abstract models of this nature were proposed before by Kosloff and coworkers 
14, lj- Here 

we consider a specific, realizable, model of a molecular level heat 



12, 



13, 



pump based on modulation of molecular energy levels. Such modulation can be achieved by 
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a stark shift affected by a tip induced .oca. e.ectric ne.d, by magnetic ne.d splitting energy 
levels and by an external force applied by the tip of an atomic force microscope |16|]. It was 
also demonstrated recently that nanotubes tension can be tuned by applying an electric field, 
thus modulating the tube vibrational frequencies |l7l.lis||. Finally, compression of molecules 
affects their vibrational modes, e.g. the radial breathing modes of nanotubes are pressure 
dependent [n| with about duo/dP ~ 1 cm _1 /Gpa [20'j, making high pressures necessary 
for a significant effect. Each of these schemes can be used as a basis of the proposed heat 
engine. Below we describe the concept of this engine, consider its performance and efficiency 
in terms of molecular and junction parameters, and suggest possible optimization methods. 



II. MODEL 



The model system consists of a molecular unit connecting two thermal reservoirs L and 
R of inverse temperatures Pl = (&b^l) _1 and [3r = (/c^Tr) -1 respectively, where k B is 
the Boltzmann constant. For simplicity we assume that heat transfer is dominated by a 
specific single mode. In addition, if the baths temperatures are low enough, only the lowest 
vibrational states of the molecule are populated, and we can model the isolated molecule 
by a two level system (TLS). An external force drives periodically the frequency of this 
molecular mode, i.e. the two level energy spacing. The total Hamiltonian therefore includes 
three terms 

H = H s + Hb + Hmb, (1) 

where 

^ = ^(|1)(1|-|0)(0|) (2) 

is the Hamiltonian of the molecular mode under consideration (h = 1). Here |0) and |1) 
represent the two states of energies e and ex, and 

LU (t) = e 1 -e = iu + F(t), (3) 

provides the time dependent driving with a static frequency ujq and a periodic modulation 
F(t) = F(t + ^y). In what follows we refer to u(t) as the instantaneous energy gap. F(t) 
can be expanded in a Fourier series 

oo 

F{t)= [A n cos(nnt) + C n sm(nttt)} . (4) 



We also define the indefinite integral of this perturbation that will be useful below 



/(f) = / F(t)dt 



'A C 

— ^ sin(nfit) cos(nfit) 

nil nil 



(5) 
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n=— oo 

The two thermal reservoirs L and R 

H B = H L + H R (6) 

do not interact directly with each other, and can exchange energy only through their coupling 
to the system. Transitions between the |0) <-> |1) states can occur due to the coupling to 
these heat baths 

Hmb = 5|0)(l| + J B t |l)(0|, 

B = B L + B R , (7) 

where B^ {K = L,R), the bath operators, are given in terms of their phonon coordinates. 
The thermal reservoirs are characterized by their spectral density functions. An essential 
ingredient of our model is having different spectral functions for the left and right reservoirs. 
Below we model this difference by assuming that the reservoirs are characterized by different 
Debye frequencies ujj ^ u^. Similar effects may be achieved by other means, e.g., connecting 
identical bathes to the system via 'doorway oscillators' of different frequencies. 

Eqs. |T|) -(jT)) represent a particular kind of a molecular relaxation process [21]. Unlike the 
standard relaxation models here the molecular mode is (a) coupled to two thermal reservoirs 
of different temperatures and spectral properties, and (b) modulated by an external force 
so that the corresponding level spacing oscillates in time. This then becomes a "driven 
dissipative" system that differs from previously considered models |22] by working in the 
system eigenstate representation and by coupling to two independent thermal reservoirs. 



III. OPERATION CYCLE 



Next we describe a setup that leads to the desired pumping operation. First, a tem- 
perature gradient is applied across the system by keeping Tr = Th > Tl = Tc- {H and 
C stand for "hot" and "cold" reservoirs). In addition, asymmetry is built into the system 
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by choosing < and k,l > kr, where Kk is a parameter related to the vibrational 
relaxation rate induced by the K thermal bath (see Eqs. (j2*Uj) and ([24))) . With this choice 
of parameters, when the TLS frequency u(t) is small, it is coupled more strongly to the left, 
cold reservoir. Energy is then injected from the left reservoir into the system whenever the 
TLS temperature defined as 

TTLs{t) = (8) 
is smaller than Tc- Here Pq and P\ are the population of the |0) and |1) states respectively. 
The TLS energy spacing is next increased by the action of the external force, therefore it 
couples more effectively to the right, hot reservoir. If the levels population is kept (almost) 
fixed during this process, the effective TLS temperature becomes very high. If it is higher 
than Tr = Tjj, heat will be transferred from the TLS into the right -hot reservoir- and the 
pumping cycle is completed. For a schematic representation see Fig. [H 

This pumping machine is a continuous version of the discrete four strokes pump of 



Ref. 



15| . Here the system is effectively disconnected from each reservoirs at different times 



due to the asymmetric construction of the reservoirs spectral properties and the system-bath 
interactions. 



IV. DYNAMICS 

Given the time dependent Hamiltonian, Eqs. (JD)-©, a Master equation for the states pop- 
ulation P n (n=0,l) can be obtained by making the following assumptions: (i) The system- 
heat bath couplings are small so that second order perturbation theory can be applied to 
yield golden-rule-type relaxation rates, (ii) The memory time of the bath fluctuations tr 
(K=L,R) is short relative to the thermal relaxation time 

r K l » T. (9) 

Here 1/r is the thermal relaxation time of the two-level system given by = [k u + kd)~ , 
see definitions below. Under these assumptions, Redfield theory leads to the Markov- 
Master equations for the states population, (See Appendix A for a detailed derivation) 

Pi = -k d (t)P! + k u (t)P ; 

Pi + Po = 1, (10) 
5 




FIG. 1: (Color online) Schematic picture of the pumping cycle, (a) At low frequencies the TLS is 
strongly coupled to the left reservoir. Thus when Ttls < Tl heat is transferred from the L bath 
to the TLS. (b) At high frequencies the TLS is coupled only to the right reservoir, thus its internal 
energy is transmitted into the right bath when Ttls > Tr. 

where 

k u = [ e iuJo{t - T) e llf ^- f{T)] (B(r)B^t))dr 
Jo 

+ f e- iuJo{t - T) e- l[m - f{T)] (B(t)B\r)}dT, (11) 
Jo 

k d = [ e iu)o{t - T) e i[f{t) -K T)] (B\t)B{T))dT 
Jo 

+ f e^^h^V®- f W(Bl(r)B(t))dT. (12) 
Jo 

Here f(t) is the time periodic function of Eq. (|5)i. Note that the relaxation rates include 
contributions from both left and right thermal baths since (B^(t)B(0)) = (B L (t)BL(0))L + 
(B' R (t)BR(0))ji, as implied by Eq. (JJJ) where the averages are over the thermal distributions 
of the corresponding baths. Therefore 

k u = k u ,L + k u , R ; k d = k d:L + k d;R , (13) 
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A. Rate constants: General expression 



We derive next explicit expressions for the rate constants in the general, non adiabatic 
regime. We begin with the first integral of the excitation rate k u in Eq. (jllj) . It can be 
expanded as follows 

/oo 
e i{uJO+nQ)x (B(0)B j (x))dx, 

(14) 

where the sum goes over Yl^ m =-oo' anc ^ ^ ne u PP er limit in the integral is extended to infinity. 
To obtain (|14|) we have utilized the Fourier expansion 



oo 

,*/(*) - 



e 

n=— oo 



J ^ mQt - (15) 



For the simple cosine modulation, F(t) = Aicos(Qt), the expansion coefficients reduce to 
the Bessel functions J n (Ai/Q) of order n. We note that the m series in Eq. (j!4j) trivially 
sums up to e % ^ l \ yet we prefer this representation since it formally eliminates the specific 
F(t) dependence from the equations. The second integral in Eq. (JTTj) can be manipulated 
in the same way to produce 

I 2 = J2 J nJ* m e i{n ~ m)nt / e i{uJo+nn)x (B(0)B\x))dx. 



n,m 



(16) 

Assuming that e i{uJ0+nU)x (B(0)B^(x))dx is a symmetric function around zero so that the 
integrals in I\ and I 2 can be replaced by | the total excitation rate becomes 

k u = h + h = Yl 38 [JnJ: n e i{n - m)Qt ] 4 n) ; 

n,m 

/oo 
e i( -" o+nn >{B(0)Bl(x))dx, (17) 
-00 

where 3? denotes the real part. It is given in terms of standard time independent transition 
rates evaluated at different overtone frequencies nu (n are integers), multiplied by the ap- 
propriate Fourier coefficients and a time dependent modulation. The downward rate kd is 
obtained in a similar way by combining the two integrals of Eq. (|12|) 

k d = J2^[Jnr m e l(n - m)m ] k d n) - 



n,m 
000 



/oo 
e i{ - U30+nQ >{B\x)B{Q))dx. (18) 
-00 



For infinitely slow modulation, A n /Q, C n /£2 — >• oo, the standard time independent expression 
for the vibrational relaxation (VR) rate |2J] is recovered, kd(t) — > by making use of the 
sum identity J2 n ,m JnJLe l(n ~ m)nt = 1- 

The transition rates (|17j) and (|18j) can be further decomposed into the L and R con- 
tributions as in Eq. (|13|) . The up and down rates induced by each thermal reservoir are 
interrelated by the detailed balance condition for each n component 

C = C^^^' (K = L,R). (19) 

Our results so far are general within the weak system- reservoir interaction limit. As a 
specific model for the bath correlation functions we invoke below the exponential energy gap 
form 25] 

, g= k. ^ +nn < u g (2o) 

I KKe -(^o+nn)/^ ; uj + ntt>uj%. 

Here is the Debye frequency characterizing the K = L,R reservoir. In general we assume 
kl 7^ kr and ^ We associate the bath relaxation time tk with the inverse Debye 
frequency. 

B. Rate constants: Adiabatic regime 

We can also derive explicit expressions for the transition rates assuming the energy mod- 
ulation is adiabatic, i.e. does not itself induce transitions in the TLS or in the thermal 
reservoirs. In this regime the integrals of Eqs. (fTT]) - (fT2*j) can be simplified by approximating 
the differences by first derivatives, f(t) — jit — x) ~ xF(t), x = t — r. Higher order terms 
are neglected assuming l^y^^l <^ (For a cosine modulation, F(t) = A\ cos(Ot), this 
condition translates into Q <C uip). Then 



oo 



h = e iuJox e i[m - f{t - x)] (B(0)B\x))dx 
Jo 

POO 

— > / e iuJ ° x e ixF{t) (B(0)B\x))dx. (21) 
Jo 

Conducting similar operations on the second integral of Eq. (fTTj) yields 

r° 

I 2 = I e iuJ0X e ixF{t) (B(Q)B^(x))dx, (22) 

J — oo 
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and the adiabatic rate constants become 

/oo 
-oo 

K,k = k^e-^^- uj(t)=u + F(t). (23) 

The adiabatic approximation therefore implies that an instantaneous detailed balance would 
be satisfied at all times if the system was coupled to a single bath, i.e. Pi(t)/P (t) = 
e -w{t)/k B T ^ w here j 1 j s the single bath temperature. Also in this case we utilize the exponential 
energy gap law for modeling the adiabatic relaxation rates 



K K , U(t) < Uj' D 

(24) 



V. ANALYSIS OF PERFORMANCE 

The performance of a heat pump can be characterized both in terms of its power- the 
amount of heat transferred per cycle (i.e. period of the modulating force) and its efficiency, 
defined as the ratio between the heat transferred and the work invested. The internal energy 
£ of the TLS is given by 

£ = e P + eiPi, (25) 

where eo and e\ are measured from some fixed reference, here chosen by Eq. (j2J) to be the 
midpoint between the two levels. The internal energy is changed either through modulation 
of the TLS energy spacing or due to population transfer between the levels [12I \\\ 

d£ dP dPi p de ^ dei 
dt dt 1 dt dt 1 dt 

This rate of energy change can be separated into its work W, and heat Q components 
Using the following equalities 

P + A=0; 60 + 4 = 0, (28) 
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that are based on Eqs. (J2J) and (fTUj). we find 

W = S(t)w(t); Q = u{t)S{t), (29) 

where S = (Pi — Po)/2 is referred to as the system polarization 12, As the effect of 
the two reservoirs is additive, we can decompose the rate at which S changes to its L and 
R contributions 

S = Sl + Sr, 

Sk = -k d , K Pi + k u , K P , (K = L,R), (30) 

where Pi and P are obtained by solving Eq. (fTU|) . Consequently, the heat flux Q can be 
written as a sum of L and R terms 

Q = Ql + Qr; Q K = u{t)S K . (31) 

We note that in steady state Ql = —Qr, i-e. the heat current is the same at the left and 
right contacts. Here these quantities are in general different, even on the average, due to 
the action of the external perturbation, J L = Ql 7^ Jr = — Qr- In the equations above the 
heat current is taken positive when flowing left to right. 

The coefficient of performance (COP) of a heat transfer machine can be defined with 
regard to its performance either as a heat pump 

rj H = Q H /W, (32) 

or as a refrigerator 

Vc = Qc/W, (33) 

where Qk = f cyde Qk {K = C, H) and W = f cyde VV. The maximal theoretical values of 
these coefficients are given by that of a reversible (Carnot) machine, 

max _ Qh _ Tjj 

Vh ~ Qh-Q c -t h -t c > (34) 

Qh — Qc Th — Tc 

In what follows we focus on the refrigerator COP, Eq. ()33|). as a measure of efficiency of our 
molecular machine. 



10 



In an ideal refrigerator the operation cycle consists four distinct steps: (i) Thermal: The 
TLS with an energy gap ojl couples to, and exchanges energy with, the left (cold) bath 
only, (ii) Adiabatic: The TLS is decoupled from the reservoirs and its energy spacing is 
increased to lur. (iii) Thermal: The TLS is coupled to, and exchanges energy with, the right 
(hot) reservoir only, (iv) Adiabatic: The TLS, again decoupled from the reservoirs, restores 
its energy gap back to the low u>l value. It should be emphasized that in the realizable 
machine discussed in Section III, system bath decoupling is not imposed. It is approximated 
by replacing the adiabatic steps by transitions whose durations are short relative to the 
thermal relaxation time associated with the system-bath coupling. Optimized performance is 
therefore obtained when the adiabatic branches of the process are fast relative to the thermal 
branches, so that no backward heat flow takes place. In contrast, the thermal branches 
should be long enough for attaining full equilibration of the TLS with the interacting bath. 

Consider first the ideal refrigerator in which the system is decoupled from the left/right 
reservoirs during the adiabatic branches. In this limit Eq. (fT0|) can be solved analytically 
for the thermal branches yielding the time evolution of the K = L, R polarization 

^(t) = ^ + (^(0)-^)e- r ^; K = L,R, (36) 

where is the equilibrium polarization 

S% = -\ tanh(u K /2k B T K ). (37) 

ujk is the time independent TLS gap when it is connected to the K reservoir and = 
kd,K + ku,K- We denote the durations of the thermal branches, i.e. the contact times of 
the system with the L and R reservoirs by tl and tr respectively, and recall that the 
polarization S = (Pi — Po)/2 does not change during the (ideal) adiabatic branches. The 
polarizations at the beginning of the thermal branches are therefore given by 5l(0) = Sr(tr) 
and Sr(0) = Sl{jl). Using these initial values in Eq. the solution of the two coupled 
linear equations at the end of the thermal branches, i.e. at t — tk is 

Sr(t r ) = S e L q + ^ t. \ (38) 

and an analogous expression for Sl(tl). 

The amount of heat pumped out of the cold (L) reservoir during each cycle is calculated 
by substituting the derivative of Sl(£) (Eq. (jHSj)) into Eq. (j2HJ), then integrating over the 
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contact time with this reservoir 



Ql = u L (S L (0)-S e L q )(e- r ^-l) 

= u L (S e L q - S%) r f~m - ■ (39) 

When the coupling times fi and f R are long relative to the inverse relaxation rates, the TLS 
equilibrates with the heat baths during the thermal branches. Then the heat pumped per 
cycle is maximized 

Q L = u L (5? - S£) . (40) 

Based on this equation we can derive the condition for attaining the desired pumping action: 
is required to be positive, implying that S e L q > S R 9 . In the classical limit, ojk < k B T K) 
using Eq. (p?7j). this translates into the condition ujl/ojr < Tl/Tr. The work performed on 
the system can be calculated similarly to yield [3| 

W = (u r -u l )(S l (t l )-S r (tr)) 

( 1 — e~ TLfL ) ( 1 — e~ rRfR ) 
= («* " -J< S ? " ^ 1 _ e -rA e -r, fa J - («) 

When TkTk — > oo, the work approaches 

W = K-^)(^-^). (42) 

The COP of this idealized machine, Eq. (JHHj) . is then 

T)c = u l /(uj r - u) L ). (43) 

It does not depend on the temperature, only on the minimal and maximal values of the 
molecular energy gap. Note that in the opposite t^Yr ~* limit, expanding e~ TKTK ~ 
1 — TrTk, leads to 

gt.Wcx /^ff* . (44) 
1 lTl + 1 RTR 

If we further assume that the contact times are proportional to the inverse of the energy 
gap modulation frequency, we conclude that both Ql and W scale like fi -1 . In the next 
section we compare these results with the performance of the realistic machine introduced 
in Section III. 
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VI. RESULTS 



In the general case Eq. (jlUj) has to be solved numerically for the populations P\(t) and 
Po(t) = 1 — Pi(t), and we use the fourth order Runge-Kutta method for this purpose. We 
focus on the long time behavior of these quantities in order to eliminate effects of the initial 
conditions. The heat current, the applied work, and the machine efficiency are calculated 
using Eqs. ()27 )) -([33 )) . In order to retain the markovian limit we choose a set of parameters 
that fulfills Up 3> T, (K = L,R). The adiabatic criteria is additionally preserved when 
du/dt < F(t)u$. 

We begin by analyzing an adiabatic machine operating under a pure sine modulation of 
the TLS gap, A n = 0, d =25 meV, C n ^ = 0, in Eq. (gj). The choice of Q=0.025 meV, 
u)p=6 meV and a;^=250 meV corresponds to the adiabatic limit. The rate constants are 
therefore calculated using equations (J2HJ)-(|2I1) instead of the general expressions (fT^) - (j2l]j) . 
This simplifies significantly the computational effort, since for C\/Q ~ 1000 expansion terms 
J n up to the order n ~1200 have to be taken into account in order to achieve convergence. We 
have also verified that the adiabatic results perfectly agree with the general formalism. The 
results of this calculation for this choice of parameters (other parameters are noted at the 
caption) are displays in Fig. |21 Shown are the TLS spacing modulation, the instantaneous 
TLS temperature and the instantaneous heat transferred at the cold and hot interfaces. We 
find that in this adiabatic limit the device does not pump heat, and energy is transferred 
from both the hot bath and the external periodic field into the cold reservoir. This is 
also demonstrated through the temperature of the TLS which is always higher than the 
temperature of the cold bath. Therefore, extraction of heat from the L reservoir is impossible, 
and Jl is negative throughout the cycle. 

Consider next the quasiadiabatic situation where the TLS energy spacing is modulated at 
a frequency that is at the same order or smaller than the inverse of the reservoirs relaxation 
times. In Fig. El we display the behavior of such a machine, where the parameters are the 
same as in the previous adiabatic model except that f2=0.5 meV. In panels (a)-(d) we show 
(a) The (periodic) time variation of the TLS energy spacing uj(t) . (b) The effective couplings 
kd,L and k^R of the TLS to the left and right reservoirs. When u(t) reaches its maximum 
value, k^L becomes negligible, ~ lxlO -4 meV, and the TLS is effectively disconnected from 
the L reservoir. In contrast, since the Debye frequency at the R side is significantly larger 
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FIG. 2: Adiabatic heat pump under a pure sine modulation, Ci=25 meV, C n ^\ = 0, A n =0. 
n = 0.025 meV, uj = 40 meV, k l = 2.5 meV, kr = 0.1 meV, lo^ = 6 meV, uj% = 250 meV, 
Tl = 200 K, Tr = 300 K. Top: Energy spacing (dashed line, left vertical axis) and the resulting 
TLS temperature calculated using Eq. (JHJ) (full line, right vertical axis). Bottom: heat current 
Jl = Ql (full) and Jr = -Qr (dashed). 

than the molecular frequencies, kd,R remains effectively constant at all times, (c) The TLS 
temperature. When u(t) becomes large, the TLS temperature reaches a maximum of ~ 600 
K, larger than Tr=300 K. Heat transfers then from the hot molecular mode to the right 
reservoir. The lowest TLS temperature of ~ 150 K is obtained at small energy spacing, 
uj=15 meV, at which k^L ~ 0.2 meV and k^R = 0.1 meV. Therefore at this point the 
molecular mode gets heat from both L and R reservoirs. This is seen in the bottom panel 
(d): Jl (full line) is positive and Jr (dashed line) is negative when the energy gap u> is in 
the neighborhood of this small value. 

For this operation mode and for these engine parameters we find that the amount of heat 
pumped out of the cold reservoir at each cycle is Q=0.23 meV, and the efficiency of the 
machine is 7/(7=0.073. Note that the heat pumped into the right contact is not the same 
due to the action of the external force. In addition, looking at the instantaneous pumping 
we observe a delay of half a cycle in the pumping action: Heat is pumped from the left cold 
reservoir when the TLS gap is minimal, and it is injected into the right hot reservoir after 
half of a cycle when u(t) becomes large. In between, due to the slow decoupling rate of the 
L reservoir from the molecule, a backward heat flow is observed. 

Figure 4 shows results pertaining to the efficiency of this machine. The top panel of Fig. 0] 
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FIG. 3: Quasi adiabatic heat pump operating under the modulation frequency £1 = 0.5 meV. Other 
parameters are as in Fig. El Shown are the TLS energy spacing modulation (a), the relaxation 
rates at the right (dashed) and left (full) contacts (b), the TLS temperature (c), the heat currents 
Jl (full) and Jr (dashed) (d). 

presents the heat transferred per cycle, plotted against the driving frequency. It reaches a 
maximal value for Q ~0.5 meV. The efficiency denned in Eq. (IHHj) increases monotonically, 
saturates, then decays very slowly. We can explain these observations as follows: For very 
slow modulation the TLS is at steady state driven by its coupling to the L and R reservoirs 
and its temperature is approximately given by T TLS ~ (X L T L + T R T R ) /(T L + T R ) [9|, higher 
than Tl = Tq. Heat then always flows towards the cold bath and Ql is negative. The 
finite coupling to both reservoirs at all times therefore inhibits the pumping operation in 
the adiabatic regime. This is in fact the extreme opposite to the optimal situation in which 
the system is decoupled from the reservoirs whenever needed, that leads to Eq. (|40p. In 
the opposite fast modulation limit the TLS temperature can reach values below T c , and 
can pump heat out of the cold reservoir as in Eq. Its efficiency is however restricted 

by the fact that time is insufficient for a full equilibration with the cold reservoir, thus the 
total energy injection is small. We have also verified that in this regime both heat and work 
decay like This implies that maximal heat pumping is obtained at some intermediate 
modulation frequency, as seen in the upper panel of Fig. 0] The machine performance is 
therefore optimal when working in the quasiadiabatic regime. 

The results obtained above make it possible to investigate ways to optimize the perfor- 
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FIG. 4: Efficiency of the non adiabatic heat pump displayed in Fig- El Top: Heat transferred per 
cycle out of the cold reservoir. Bottom: The refrigerator coefficient of performance (Eq. (|33j0 . 
Ci =20 meV (dashed), C x =25 meV (Full), d =30 meV (dashed-dotted). 
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FIG. 5: (Color online) A surface plot of the heat pumped displayed against the driving frequency 
Q and amplitude C\ for a system characterized by the parameters of Fig. 01 

mance of our heat pump. Figure El presents a surface plot of Ql per cycle as a function 
of amplitude C\ and frequency Q for a sine type perturbation. We find that for very weak 
modulations the system cannot pump heat, i.e. Ql is negative. The maximal amount of 
heat per cycle of Ql=0.4 meV is pumped at C\ ~ 30 meV and fl ~ 0.5 meV. 

Another technique for optimizing the heat pump operation is by devising an optimized 
shape for the modulation function F(t). As discussed above, F(t) should be designed so as 
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FIG. 6: Non-adiabatic heat pump. C n = 4/Vn x 18 meV, n=1..9, A n =Q, = 0.5 meV, ujq = 40 
meV, k l =2.5 meV, k r = 0.1 meV, uj^ = 6 meV, = 250 meV, T L = 200 K, T R = 300 K. (Top) 
Shown are the TLS energy spacing modulation (left) and the resulting TLS temperature (right). 
The bottom panel displays the heat currents Ql (full) and —Qr (dashed). 



to minimize reverse heat transport processes. The parameters that can be manipulated are 
the functional form of the modulation, the total time duration of the pulse, and the time 

n 

allocated to the four operation branches |15| . 

As displayed in Fig. 01 when utilizing a sine modulation with the given amplitude and 
coupling parameters, the TLS temperature varies between the minimal value of T=150 K 
< Tl and the maximal value T=600 K > Tr. We have argued that the machine efficiency 
can be improved if the thermal branches are long enough so that the TLS comes as close 
as possible to equilibrium with the corresponding thermal bath, and if between the thermal 
branches the energy gap is varied as rapidly as possible. In Figure |H1 we show such a machine 
where the modulation is tailored such as to "wait" at its maximal and minimal values giving 
the TLS more time to equilibrate with the different reservoirs at different parts of its cycle. 
At the same time the energy gap is changed rapidly in order to eliminate backward flow 
to the L bath. The top panel (dashed line) presents the shape of the driving signal and 
the ensuing TLS temperature. This modulation function is constructed from the series 
F(t) = ^2 n=1 3 L C n sin(fint), C n — 4/irn x 18 meV, L n =9. We find that indeed during the 
"wait" time at the minimal energy spacing ~ uj = 20 meV the TLS heats from 150 K up to 
~200 K. The same effect is observed when the TLS reaches its maximal u = 60 meV value: 
The TLS temperature reduces from ~ 600 K down to ~ 400 K. The fast oscillations of Jl at 
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high temperatures (e.g. between 40-45 ps) sum up close to zero. The refrigerator efficiency 
is 770=0.021 and <2l=0.06 meV per cycle. The corresponding machine operating under a 
pure sine modulation with the same amplitude (L n =l, C%=18 meV) does not pump heat. 

VII. SUMMARY 

In analogy with electron and spin pumps that were investigated in recent years, we 
propose here a molecular level thermal machine that pumps energy from a cold to a hot 
reservoir. We discuss its operating principles, performance, optimization techniques, and 
possible physical realizations. We expect that in the future devices whose functionality is 
determined by both their thermal and electrical properties j3] will be of great interest. 

APPENDIX A: Derivation of the quantum master equation. 

Here we derive the quantum master equation for a two level system with a time dependent 
energy spacing which is interacting with two thermal reservoirs, Eqs. (HJ-Q 



H B = H L + H R , and the system bath coupling B includes the L and R terms, B = B L + B R . 
The evolution of the total density matrix is given by the Liouville equation (h = 1) 



The equations of motion for each density matrix component pij are given in terms of the 
bath operators 



H = ^o + F(t))(|l>(l|-|0><0|) 
+ 73|0)(l| + 73 t |l)(0| + ifB. 



(Al) 



dp 
dt 



i[H,p]. 



(A2) 



Pi,o 



Po,o 



Po,i 



-iB*po tl + ipiflB, 

-iBp lj0 + ip 0:1 B f , 

i (u + F(t)) p ,i - iBpx,i + ipo,oB, 

-i (uq + F(t)) px,o - iB^p 0fi + ipi tl 



St. 



(A3) 
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Next we formally integrate the nondiagonal terms po,i an d pi,o using the Leibnitz integral 
rule 



dt 



d r v ® 

- I fit, r)dr = v'(t)f(t, v(t)) - u'(t)f(t, u(t)) 

u{t) 

Mt) o 

+ J t g- t f(t,r)dr, (A4) 



and obtain 



Po,i(f) = J e- ( t -^e i ^W-^^[-^(T)p lil (T)+zpo,o(T)5(r)]dT. 

Pi,o(t) = pg,i(0 (A5) 



where f(t) = f F(t)dt. We substitute these expressions into the equations of the diagonal 
terms po,o and and trace over both L and R thermal baths assuming the density matrix 
can be decomposed at all times by p(t) = PLPR&if)- Here a is the reduced density matrix 
operator and px = e~" KHK /Trace(e _/3Ki ^ K ), K = L,R. Following the standard Redfield- 



Bloch derivation 23|, i.e. second order perturbation theory combined with the assumption 
that bath correlation functions decay rapidly on the time scale of the change of a, we obtain 
the quantum master equation for the diagonal reduced density matrix elements P n = cr„ in , 

71=0,1 

ft 



p 1 = P (t) / e iuJo{t - T) e i[m - f{T)] (B(r)B\t))d 



r 



+ P (t) / e-^-^e-^^^iB^B^r^dr 



o 



Pi(t) / e^ o{t - T) e i[m ^ f{T)] (B\t)B(T))dT 



- P^t) / e^^e^V®-tW(B\T)B(t))dT\ 
Jo 

P (t) = l-Pi(t). (A6) 

In the markovian limit we further extend the upper limit in the integrals to infinity, and 
assume that bath correlation functions do not depend on the initial time. Note that no 
restrictions are imposed on the modulation term F(t), e.g. in general it need not be periodic. 
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